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The coarsening dynamics of the Cahn-Hilliard equation with order-parameter dependent mobility, 
\{(j>) oc (1 — 4> 2 ) a , is addressed at zero temperature in the Lifshitz-Slyozov limit where the minority 
phase occupies a vanishingly small volume fraction. Despite the absence of bulk diffusion for a > 0, 
the mean domain size is found to grow as (R) oc t 1 ^ z+a \ due to subdiffusive transport of the order 
parameter through the majority phase. The domain-size distribution is determined explicitly for 
the physically relevant case a = 1. 
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The phenomenon of spinodal decomposition in, e.g. 
binary alloys, is usually modelled by the Cahn-Hilliard 
equation (!]] 

£-V.(AV,)-V.(»v£), (D 

for the order-parameter field <fi. Eq. ([j]) takes the form 
of a continuity equation, dt<j) = — V • j, with current j = 
— AV/x, where A is a transport coefficient ('mobility') and 
the chemical potential \i is the functional derivative, /x = 
5F/5<fi, of a Ginzburg-Landau free energy functional F[tfi] 
given by 

F[cf>} = J d d x(±(V^ + V(<j>)) . (2) 

Here V((f>) is the usual double-well potential whose min- 
ima (taken here to be at <f> = ±1) represent the equilib- 
rium phases. 

In conventional treatments of ([!]), the mobility A is 
taken to be a constant, i.e. independent of the order pa- 
rameter <j). Recently, however, there has been consider- 
able interest in cases where A depends explicitly on 
(j), notably through the dependence A(0) = Ao(l — </> 2 ). 
This interest has a physical origin. It has been noticed ||] 
that when one models the coupling to an external driv- 
ing field E, such as gravity, through an additional term 
Fx[(f>} = -E J d d xz4>(x) in F[(j>] (where here the field E 
is in the z-direction) , this extra term does not change (|l|) 
unless A depends on cA. This is because 5Fi/5(f> = —Ez, 
and V 2 z — 0. Physically, it is clear that an external field 
of this form accelerates the phase separation, so A must 
be 0-dependent. Indeed, phenomenological derivations 
P,pl of A yield precisely the form A oc (1 — <fi 2 ) alluded 
to above. Furthermore, the coarsening dynamics of this 
model has been studied by computer simulations, both 
with (fH^II and without || external driving forces. It 
is therefore interesting to study this problem in its own 
right, even without external driving forces. 



In this Communication, we study a general class of 
systems described by Eq. ([j]) with 

A = (i-^ 2 r (3) 

(we absorb the constant Ao into the timescale). To make 
analytical progress, we specialize to the case where the 
minority phase occupies a vanishingly small volume frac- 
tion. For the conventional case (a = 0), this is the 
limit treated by the seminal work of Lifshitz and Sly- 
ozov (LS) j|o|, and by Wagner pH , which leads to the 
result (R) oc t 1 / 3 for the mean domain size, and gives 
an exact expression for the domain-size distribution. For 
general a > we find (R) oc i 1 /( 3+Q ). We also determine 
explicitly the domain-size distribution for the physically 
relevant case a = 1. (The other physically relevant case, 
a = 0, has been treated by LS.) 

For small volume fractions, coarsening proceeds by nu- 
cleation and growth rather than by spinodal decomposi- 
tion. For present purposes we limit discussion to the 
late stages of growth, described by the LS evaporation- 
condensation mechanism, by which large domains (of ra- 
dius R > R c (t)) grow at the expense of small ones (with 
R < R c ). In this regime the droplet size distribution has 
the scaling form n(R) = R~ 4 f(R/R c ). 

Before we proceed, two comments are in order. We 
have noted that phenomenological considerations indi- 
cate that (|^) with a = 1 is an appropriate form for the 
mobility in the Cahn-Hilliard equation ([!]). This being 
the case, one may wonder about the physical relevance 
of the conventional (i.e. with a = 0) Cahn-Hilliard equa- 
tion. The point is that thermal noise, omitted from (Q), 
reduces the magnitude of equilibrium order parameter 
from unity to (f>Q. Since thermal fluctuations are irrel- 
evant on large scales |l]], however, one can continue to 
work with the noise- free equation ([!]), provided that (in- 
ter alia) one uses a renormalized potential whose min- 
ima are now at ±4>o- The bulk mobility then takes the 
value Xbuik = 1 — 4>o > 0, and conventional LS behavior 
is recovered. The relevance of a = 1 is then limited to 
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'deep quenches', where the effect of thermal noise is small 
enough that the predicted t 1 ! 4 growth extends over a sig- 
nificant time domain (before t 1 / 3 LS growth eventually 
sets in). In simulations, of course, one can simply work 
at zero temperature, when the t 1 / 4 behavior (or iVw+oO 
in general) will describe the asymptotic growth. 

The second comment concerns the role of surface dif- 
fusion. It has often been stated that (without thermal 
noise) a factor (1 — (j> 2 ) in the mobility prevents bulk 
diffusion, and therefore surface diffusion (i.e. diffusion 
along the interfaces), which leads to t 1 / 4 growth, is the 
dominant coarsening mechanism in this case. It is true 
that the bulk diffusion constant vanishes for a > [see 
@ below]. In the far off-critical systems discussed here, 
however, where the minority phase does not percolate, 
surface diffusion alone cannot lead to large-scale coars- 
ening. It turns out that for a > there is still bulk 
transport, although this of a subdiffusive, rather than dif- 
fusive, character. 

We begin by considering a single spherical domain of 
'plus' phase, with radius R, immersed in a sea of 'minus' 
phase. We suppose that the minus phase is supersatu- 
rated with the plus phase, i.e. <fi = — 1 + e at infinity, 
and we work in the limit of small supersaturation, e <C 1 • 
First note that the chemical potential /i is given by 
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(4) 



where the prime indicates a derivative. In the bulk 
phases, away from the interface, <p varies slowly in space 
and the V 2 </> term in (g) can be neglected. Setting 
cj> = — 1 + <j> in the minus phase (with cj> = e at infin- 
ity), Eq. (Eh gives, to lowest order in </>, 



/*(r) = V"(-l)4>(r) 



r>R 



(5) 



Inserting this result into ([j]), using (||) for A, and retaining 
the leading terms for (j> <C 1, gives, in the minus phase 
(away from the interface), 



d t <t> = 2 a V"(-l)V • (<TV0) . 



(0) 



Note that, except for a — 0, this equation is not of the 
usual diffusive form. We shall find, nevertheless, that 
it still leads to bulk transport, albeit of a subdiffusive 
character. 

We now make the usual assumption (to be verified a 
posteriori) that the interface moves slowly enough (for 
large R) that the chemical potential is always in equilib- 
rium with the interface. Then the time-derivative term 
can be set to zero in (||). Using the linear relation 
between /j, and <j) in the bulk minus phase, (g) can be 
recast as 



V 2 (fi i+a ) = 



(7) 



a simple generalization of the Laplace equation V 2 /i = 
obtained when a = 0. 



What are the boundary conditions on (g)? At infinity, 
we have 



Ai(oo) = A^oo - V"(-l)£(oo) = V"(-l)e 



(8) 



The second boundary condition, at r = R, is just the 
usual Gibbs-Thomson boundary condition 



fi(R) = a/R , 



(9) 



where a is the surface tension. To derive (g), one 
first writes (exploiting the spherical symmetry) V 2 = 
d 2 (f> + {2/r)d r (f> in (g). Then one multiplies (g) by d r <p 
and integrates across the interface. Using the fact that 
d r 4> is sharply peaked at the interface gives n(R)A<j) = 
-(2/R) J dr(d r 4>) 2 = -2a /R, where A0 = -2 is the 
discontinuity in across the interface. This reproduces 

The solution of (m) with boundary conditions (|3|) and 
(I) is 



a\ 1+a 

r) 



R 



r > R 



(10) 



The time-dependence of R is obtained by considering 
the flux of material to (or from) infinity. The current 
j through the minus phase is 



j(r) = -Xd r n= -(24>) a d rf i 



2" 



{l + a)[V"(-l)} a 



/ a\ 1 + a 

\r) 



R 



(11) 



leading to an outward flux of material / = — 2 x 
4nR 2 dR/dt = ATrr 2 j(r) (where the factor of 2 on the left 
represents the difference of between the two phases). 
This gives 



dR C 



dt 



R \R. 
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where C = [2cr/U"(-l)] a [cr/2(l + a)} is a constant, and 
R c = a I '/ioo = a/eV"(— 1) is the critical radius, i.e. the 
domain will grow if R > R c and shrink if R < R c . 

For the case of zero supersaturation (e = 0), R c = oo 
and dR/dt = —C/R 2+a . In this case all drops shrink 
(by evaporation of material to infinity). The collapse 
time t c for a drop of initial size R scales as t c cx R 3+a , 
which already suggests the scaling R c (t) ~ t 1 /( 3 + Q ) when 
evaporation and condensation mechanisms compete in 
the many-domain situation. 

Consider now a dilute assembly a spherical drops of 
various sizes. The derivation, for general a, of the scal- 
ing distribution of sizes follows that of LS for a = 0. 
The basic idea is that one has a time-dependent critical 
radius R c (t) which is determined self-consistently. Sup- 
pose that, in the late stages of growth, the distribution 
of domain radii is given by the scaling form 
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n(R,t) 



I 
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where n(R,t)dR is the number of domains per unit vol- 
ume with radii in the interval [R, R + dR). The pref- 
actor R~ ensures that the volume fraction occupied by 
the domains, ip — J dR (4ttR 3 /3)n(R, t), is conserved. 
Inserting (O) into the continuity equation 



dn d 



+ — («n) = 0, 



dt dR 



(14) 



where v — dR/dt is the domain- wall velocity given by 
@, yields 



c 



R. 
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(4/ + xf) 



f' + 



„3 + Q 
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where R c = dR c /dt. Consistency requires that the R c 
dependence drop out of this equation. This gives 



Ri +a R c = >yC , 
where 7 is a constant, implying 

R c (t) = [(3 + a) 7 Ct] 1/(3+Q) 



(16) 
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Thus the characteristic domain size grows as t 1 /( 3+Q ) as 
anticipated. This result generalizes the usual t 1 ^ 3 LS 
growth law. 

Using ( |l7|) in ( |l5|) , the latter can be integrated directly, 
in the form 



( 1 
f 



dx 2 

x 



a — x 



l+a 



47X 



3+Q 



JX 



3+a ™l+a 
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(18) 



where we remind the reader that x is the scaled radius, 
x = R/R c . 

It would seem that there is a family of solutions, pa- 
rameterized by 7. In fact this is not so - there is a unique 
value of gamma, determined following the method used 
by LS |T(| . First we recall that conservation of the order 
parameter requires that the total volume of the domains 
in the late-stage scaling regime (where the value of the 
order parameter in the majority phase approaches —1) 
be conserved, i.e. 



A roo a />oo 

— / dRR 3 n{R,t) = — dxx 3 f(x)=ip, 
3 Jo 3 Jo 



(19) 



where tjj is the volume fraction of the minority phase. It 
follows that there is a maximum value, x m , of x above 
which / must vanish. Otherwise, ( |Ts| ) implies / ~ x~^ 
for x — ► 00, and the integral ( |l9| ) will be (logarithmically) 
divergent. In fact, the denominator of the integral in ( |l8| ) 
must have a double zero at x = x m . To see why this is 



so, consider the time evolution, for a given domain, of 
the scaled radius x = R/R c . From dT2n one obtains 



R RR C 



R c Rl 
l fl 
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(3 + a)jb \x x 2+a 

(3T^ 5(x) ' 
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(20) 
(21) 
(22) 



the last equation defining g(x). 

The function g(x) has a single maximum in the interval 
(0, 00). If 7 < 70, where 
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l + a 



3 + a J \3 + a: 



2/ (l+a) 



(23) 



then the maximum lies above the x-axis and g(x) has two 
zeros x\ and X2, with xi < X2- Under the dynamics (£2), 
if a; < xi initially, then x flows to zero, whereas if x > x\ 
then x flows to x%. However, as t — > 00, X2R C —> 00 
which violates the conservation of the order parameter. 
Similarly, if 7 > 70, g(x) is negative everywhere in (0, 00), 
and all domains flow to zero size, which again violates the 
conservation. We conclude that 7 = 70 

For 7 = 70, the maximum of g(x) lies on the x axis, 
and g{x) has a double zero at 



1/(1+") 



Eq. (g) then gives 

2 + a - x a+1 - 4j Q x 3+a 



lnf(x) 



x(j x 



3 + Q 



1) 



dx 



f(x) = , otherwise. 



(24) 



x < x m , (25) 
(26) 



While this integral cannot be evaluated for general a, the 
scaling function can be determined for the two cases of 
greatest physical interest. For a = 0, 

f(x) = aoi'x 2 (3 - 2xy 11/3 (x + 3)" 7/3 exp (- 



3 - 2x 
(27) 

for x < x m — 3/2, and the usual LS result |l(| is re- 
covered. The normalization condition (19) gives ao = 
186.13.... 
For a = 1. 



f(x) = a^x 3 (2-x 2 )" 7/2 exp ^~ 



(28) 



for x < where ai = 7.785 .... Eq. ( |2^ ) gives the scal- 
ing distribution of domain sizes for the modified Cahn- 
Hilliard equation simulated in reference ||, although in 
that work the domain-size distribution was not measured. 



3 



It remains to justify the claim that one need consider 
only stationary solutions of (6), i.e. that the interfaces 
move so slowly that they can be regarded as station- 
ary while cf> relaxes. In other words, we want to justify 
the 'adiabatic' approximation of treating (f> as given by 
its equilibrium configuration for the instantaneous posi- 
tions of the interfaces. We consider late-stage coarsening, 
when R c is the only characteristic length scale. Using cf> = 
fi/V"(-l) ~ a/R c V"(-l), Eq. (6) gives the relaxation 
time of cf>, for fixed interfaces, as t re i ~ R^/V" (— l)(f> a ~ 
R 2 + a V"(-l) a - 1 /a a . From (|l|), the characteristic in- 
terface velocity is dtR c ~ C/R 2+a , where C is given 
below Eq. (|lj) . This gives the typical distance moved by 
an interface in time t re i as t re i &tR c ~ o~/V"(— 1). This 
fixed length (of order the interface thickness) is negligible 
compared to R Cl justifying the approximation. 

We note that the derivation of the i 1 /( 3+Q ) growth re- 
quires that the potential V{(j>) have quadratic minima, 
since the amplitude of the power-law growth depends 
(except for a = 0) on V"(—l) through the constant C in 
(|l7|). It is interesting that for a — 0, R c {t) depends on 
the potential only through the surface tension a, and is 
therefore independent of the detailed form of V(4>). The 
question of the appropriate form of the potential for deep 
quenches deserves further consideration. 

We stress that the results presented here are, like the 
original LS calculation, valid only in the limit where the 
minority phase occupies an infinitesimal volume fraction. 
In particular, we anticipate significant corrections to the 
domain size scaling function ( p8| ) even for ip as small as 
1CP 2 . This is certainly so for a = 0, and improved forms 
for f(x) have been suggested by a number of authors (L2] . 
It would be interesting to see whether similar techniques 
can be used for general a. 

By contrast, we expect the domain size growth law, 
(R) ~ i 1 /( 3 + Q ) to hold whenever the minority phase 
consists of isolated domains. Lacasta et al. jl| measure 
an effective growth exponent of 0.20 ± 0.01 from two- 
dimensional simulations with a = 1 and ip — 0.3. The 
extrapolation to late times (their Figure 3) required for 
comparison with the present predictions does not, how- 
ever, seem completely straightforward. Also it should 
be noted that the present theory is restricted to three 
dimensions. While extension to general d > 2 is straight- 
forward, and will not change the growth law, the case 
d = 2 is special (because of the singular form of the 
Green's function for the Laplacian) and requires a sepa- 
rate study |0 . 

If both phases are continuous, one needs to consider the 
competing surface diffusion process, which leads to t 1 / 4 
growth for the characteristic length scale (defined now, 
for example, by the first zero of the pair correlation func- 
tion). Wc conclude that, when both phases percolate, 
bulk transport (i 1 /( 3+Q ) growth) dominates for a < 1, 
surface diffusion (i 1 / 4 growth) dominates for a > 1, while 
both processes contribute the same i 1//4 growth for a = 1. 



For ip = 0.5 (and a — 1) Lacasta et al. Q measure an 
effective growth exponent of 0.22 ± 0.01, again below the 
expected value of 1/4. Longer runs may be helpful in 
clarifying whether this discrepancy is real. 

In summary, the Lifshitz-Slyozov theory of late-stage 
coarsening has been generalized to a class of models with 
vanishing bulk mobility. Coarsening occurs by subdiffu- 
sive transport of the order parameter through the major- 
ity phase. The growth law for the mean domain size, and 
the scaling form for the domain-size distribution, have 
been determined in the LS limit where the minority phase 
occupies a vanishingly small volume fraction. The result 
for the growth law, however, should hold whenever the 
minority phase consists of isolated domains. 
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